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ABSTRACT 

The five millisecond pulsars that inhabit NGC 6752 display locations or accelerations that are quite 
unusual compared to all other pulsars known in globular clusters. In particular PSR-A, a binary pulsar, 
lives in the cluster halo, while PSR-B and PSR-E, located in the core, show remarkably high negative 
spin derivatives. This is suggestive that some uncommon dynamical process is at play in the cluster 
core that we attribute to the presence of a massive perturber. We here investigate whether a single 
intermediate-mass black hole, lying on the extrapolation of the mass A^bh versus a relation observed 
in galaxy spheroids, or a less massive binary black hole could play the requested role. To this purpose 
we simulated binary-binary encounters involving PSR-A, its companion star, and the black hole(s). A 
"stellar-mass binary" black hole of (5OM , 10M Q ) can imprint the right thrust to propel PSR-A in the 
halo during a flyby. The flyby is gentle and does not alter the internal properties of the binary pulsar. An 
"intermediate-mass binary" black hole of (200M Q , 10M Q ) tends to impart a recoil speed larger than the 
escape speed: It can release PSR-A on the right orbit if its separation is wide. A "single intermediate- 
mass" black hole of mass A^bh ^ 500M Q may have ejected PSR-A at the periphery of NGC 6752 in a 
close dynamical encounter involving the binary pulsar, the black hole and a star belonging to its cusp. 
The encounter gives correct speeds but alters significantly the eccentricity of the impinging binary, so 
that it must occur before the neutron star of PSR-A is recycled via accretion torques. The influence of 
an intermediate-mass binary black hole on the acceleration of the two core pulsars is studied, and the 
ejection of stars by the binary is shortly explored. In inspecting our close 4-body encounters, we have 
found that a single or binary black hole may attract on a long-term stable orbit a millisecond pulsar. 
Timing measurements on the captured "satellite" pulsar, either member of a hierarchical triple or of the 
cusp, could unambiguously unveil the presence of a black hole(s) in the core of a globular cluster. 

Subject headings: stars: neutron - stars: black hole - pulsars: general - stars: globular clusters 

1. INTRODUCTION 

Five millisecond pulsars have been recently discovered in 
NGC 6752 displaying unexpected characteristics (DAmico 
et al. 2002). PSR-A, a canonical recycled binary pulsar, 
holds the record of being the farthest millisecond pulsar 
ever observed from the gravitational center of a globular 
cluster, at a distance of ~ 3.3 half mass radii. PSR-C, an 
isolated pulsar, ranks second in the list of the most offset 
pulsars, being at a distance of 1.4 half mass radii from the 
center 1 . PSR-B, PSR-E, and PSR-D are located instead 
within the cluster core and are single; PSR-B and PSR-E 
have remarkably high negative spin derivatives while PSR- 
D has a positive P, one of the highest measured among the 
globular cluster pulsars. If the negative spin derivatives of 
PSR-B and PSR-E are ascribed to the overall effect of the 
cluster gravitational potential, this would result in a cen- 
tral projected mass to light ratio M/Ly ^6 — 7 (Ferraro 
et al. 2003a). This high ratio may imply the presence of 
« 1,000 — 2,000M Q of under- luminous matter enclosed 
within the central 0.08 pc of the cluster (Ferraro et al. 
2003a), perhaps in the form of a single intermediate-mass 

1 D'Amico et al. (2002) considered the possibility of a chance superposition of a galactic field millisecond pulsar inside the observed projected 
distance from the center of NGC6752, inferring a probability < 10 — 4 for each of the two offset sources. 
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black hole and/or in the form of collapsed stellar remnants. 
In this scenario even the very high positive P of PSR-D 
could be explained as due to the line-of-sight gravitational 
pull of this unseen matter. However, as pulsars spin down 
due to rotational energy losses, the period derivative of 
PSR-D could also be intrinsic. Alternatively, one could 
argue that the negative values of P of PSR-B and PSR- 
E would result from the gravitational pull of some local 
perturber, such as a nearby passing star or even a more 
massive exotic objects (Ferraro et al. 2003a). The millisec- 
ond pulsars in NGC 6752 are peculiar in their location or 
acceleration, and the combination of these facts strongly 
suggests the occurrence of an uncommon dynamics in the 
core and halo of NGC 6752. It is our aim to address this 
issue here in detail. 

In a previous paper (Colpi, Possenti, & Gualandris 2002, 
CPG hereafter) we explored a number of roots for the ori- 
gin of PSR-A: PSR-A may have originated from a primor- 
dial binary, born either in the halo or in the core; but a 
careful analysis (based on considerations on characteristic 
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lifetimes and neutron star natal kicks) led us to discard 
both these two hypotheses (see CPG). A third, involving 
a 3-body scattering or exchange event off core stars, was 
also rejected given the tight constraints imposed by the 
binary nature of PSR-A (CPG). We thus were led to con- 
jecture that a more massive target, such as a binary of two 
black holes with masses in the range w 10 — 100M Q could 
have provided, in a 4-body scattering event, the necessary 
thrust to propel PSR-A into its current halo orbit, at an 
acceptable event rate (CPG). 

PSR-A may just signal the presence of a black hole bi- 
nary in NGC 6752, but this is somewhat puzzling. Black 
hole binaries are expected to form in rich star clusters; 
but along the course of evolution they are expected to exit 
their parent cluster (Kulkarni, Hut, & McMillan 1993; Sig- 
urdsson & Hernquist 1993). The black holes (relics of the 
most massive stars) tend to pair with other black holes in 
binaries as soon as they segregate in the cluster core by dy- 
namical friction (Kulkarni, Hut, & McMillan 1993). Some 
binary may rapidly merge emitting gravitational waves 
(Benacquista 1999; Miller 2002) creating a more massive 
black hole (Miller & Hamilton 2002), but most/many leave 
the cluster since close (3- or 4-body) dynamical encounters 
among the black holes eject them, single or in binaries, due 
to recoil (Sigurdsson & Hernquist 1993; Kulkarni, Hut, & 
McMillan 1993). Current N-body simulations (Portegies 
Zwart & McMillan 2000) expect no black hole or one black 
hole binary to remain in the cluster. Thus, if the interpre- 
tation of PSR-A is correct, we may provide first dynamical 
evidence of the "only" black hole binary that avoided the 
escape. A goal of this paper is to further explore this possi- 
bility narrowing the range of masses of the two black holes 
from pulsar dynamics. 

There is an alternative root for the origin of the un- 
usual location of PSR-A that we wish to explore here in 
connection to the high M/Ly ratio, and so to the potential 
presence of under-luminous matter in the cluster core: the 
ejection of PSR-A into the halo of NGC 6752 by a dynami- 
cal encounter with a central intermediate-mass black hole. 
According to an old suggestion by Frank & Rees (1976) 
stars in the cusp of a massive black hole can eject other 
stars plunging inside (Lin & Tremaine 1980). A flyby in- 
volving a central black hole, a bound star orbiting around 
it, and the binary pulsar, may imprint to PSR-A the right 
recoil speed to climb the potential well of the cluster and 
reach the halo. 

Single black holes of intermediate mass (A4bh^5OOM0) 
may inhabit the center of globular clusters; their formation 
root can involve the runaway growth of a super-massive 
star through collisions of heavy stars in the young cluster 
(Portegies Zwart & McMillan 2002) or the occurrence of 
repeated mergers among compact objects (Miller & Hamil- 
ton 2002; see van der Marel 2003 and Miller 2003 for a 
review). Recently, HST/STIS observations of the glob- 
ular cluster Gl in M31 (Gebhardt et al. 2002) and of 
M15 in the Milky Way (Gerssen et al. 2002, 2003) have 
provided first clues for the presence of a single central 
black hole. Surprisingly, the two postulated black holes, 
of mass 2.0:^8 x 1O 4 M in Gl, and of 1.7± 2 { 7 7 x 1O 3 M 
in M15, seem to lie just along the extrapolation of the 
black hole mass versus ID dispersion velocity relation 
(A4bh — f) obeyed by the super-massive black holes in 



galaxy spheroids (Ferrarese & Merritt 2000; Gebhardt et 
al. 2000; Tremaine et al. 2002). Here we assume that 
the hypothetical black hole in NGC 6752 just lies on the 
A4bh — & relation. (Note that current evidence of a large 
black hole in Gl and M15 is not yet compelling: Baum- 
gardt et al. 2003a, and 2003b provide an alternative to 
the black hole hypothesis, i.e., a cluster of lower-mass col- 
lapsed objects.) 

In this paper, we simulate a series of binary-binary en- 
counters to study (i) the scattering of the binary pulsar 
PSR-A off a stellar mass or intermediate-mass black hole 
binary, and (ii) the plunge-in of PSR-A inside the cusp 
of a single intermediate-mass black hole. In §2 we select 
the mass spectrum of the two black holes in the binary and 
choose the mass of the single intermediate-mass black hole. 
In §3 we describe properties of the binaries, and identify 
the possible end-states of binary-binary encounters (flyby, 
ionization and exchange). The binary hosting the neu- 
tron star is described in its post-recycling phase as well as 
in its pre-recycling state. In §4 we extract from a series 
of binary-binary simulations post-encounter distributions 
of recoil velocities, eccentricities and binary separations, 
running more than 20, 000 encounters. In §5 we examine 
the statistics of the encounters, while in §6 we shortly ex- 
plore exotic end-states. In §7, we model the random walk 
(moderated by dynamical friction) of the black hole binary 
induced by stellar collisions, and try to evaluate the local 
gravitational pull that a black hole binary may exert on 
the core pulsars, inside NGC 6752. In §8 we explore the 
survival of a black hole binary in NGC 6752, mimicking 
mass growth and hardening by stellar encounters. In §9, 
we outline the key findings and present our conclusions. 

2. BLACK HOLE MASSES 

The masses of the two black holes in the binary, de- 
noted Mbh (for the heavier) and tobh (for the lighter), 
are a free parameter (let Mt = Mbh + wbh be the total 
mass of the binary). Constraints on their values derive 
mainly from limits on their ejection due to gravitational 
encounters and survival against coalescence by emission of 
gravitational waves (GWs hereafter). 

In general, survival inside the cluster (followed by coa- 
lescence due to GWs) would win over ejection if the bi- 
nary, at the separation a eq where the time-scale for coa- 
lescence t gw = [5c 5 a4 q (l - e 2 ) 7 / 2 ]/[256G 3 MB H ™BHM T ] 
becomes comparable to the collision time tbh ~ 
0-BH/[nBHGM T 7ra O q], has a binding energy E hin = 
GMBH"T.BH/2a oq smaller than the minimum binding en- 
ergy for expulsion, _E j(V^ s ), given an encounter with a 
field black hole (here, «bh and <7bh are the density and 
ID dispersion velocity of the field black holes interacting 
with the binary, and V cs is the escape speed from the core) . 
If rriBH,f is the mass of the field black hole impinging on 
the binary, momentum conservation implies a minimum 
energy for escape equal to 

E cj (V cs ) = (-M (M T + m BH ,f) (—) V e \, (1) 

\ 2 ?BH/ \ m BH,f/ 

where £bh controls the relative energy exchange per scat- 
tering (see Table 1 and Table 3 for the values of £bh)- If 
-E c j > -Ebin at a oq , the binary avoids ejection during its en- 
tire lifetime (when its separation a > a cq ), and may even- 
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tually merge producing a more massive black hole when a 
drops below a cq . 
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Fig. 1. — A^BH.crit against meH (in units of Mq), for a close 
dynamical encounter with a field black hole of mass ttibh f = fBH 
(black strip) and mgn f = IOMq (dark grey strip). The binary 
has eccentricity e = 0.7. The escape speed is V cs = 40kms — 1 
and £bh = 0.5 (Miller & Hamilton 2002). The background black 
hole density, and dispersion velocity are set equal to 10 6 pc -3 and 
7kms — 1 , respectively. In generating the strips we introduced some 
scatter in §bh> V e s and density, to bracket uncertainties. The shad- 
owed areas in light grey above each strip indicate the permitted 
values of AfBH> consistent with binary black hole retention in the 
cluster. 

The condition E C j(V cs ) = Eb in (a cq ) selects a critical 
mass AfBH.crit for the heavier black hole in the binary, for 
fixed tobhj above which ejection by recoil off a field black 
hole is avoided. In Figure 1, we plot A/BH,crit against tobh, 
considering a close dynamical encounter with a field black 
hole of mass TOBH,f • The black strip refers to a mass msH.f 
set equal to tobh, while the dark grey strip refers to scat- 
tering with a field black hole of fixed mass m.BH,f = WMq. 
Above the two strips (shadowed light grey areas depend- 
ing on JTiBH.f ) binaries with Mbh > MeH.crit remain in the 
cluster. 

Figure 1 shows clearly that black holes with larger 
masses (>30M Q ) and small mass ratios tobh/Mbh are 
preferentially retained. We have considered values of 
(M BH /Mq, to bh /Mq) equal to (10,10), (30,3), (50,10), and 
(200,10) for the binary black holes in our simulations of 
binary-binary encounters; in one series of run we consid- 
ered a (50,1.4) binary comprising a black hole and a canon- 
ical neutron star 2 . The (10,10) binary is in the forbidden 
region; it is studied to set comparisons with heavier bina- 
ries, and because its mass is in the range of observed black 
hole masses (Bailyn et al. 1998). Binaries with the heavi- 
est black hole below 100M Q will be considered as stellar- 
mass black hole binaries. Those with Mbh > 100M Q will 
be referred to as intermediate-mass black hole binaries. 




1000 



100 



50 100 

cr (km s _1 ) 



500 



Fig. 2. — .Mbh versus a relation for super-massive black holes 
in galaxy bulges (Tremainc et al. 2002). Filled dots refer to the 
observation for Gl, and M15, and to the hypothetical black hole in 
NGC 6752. The horizontal error bar for the black hole in NGC 6752 
indicates the uncertainties in the measure of the central stellar dis- 
persion velocity. The upper horizontal line refers to the upper limit 
on under-luminous matter present in the core as inferred from the 
pulsar acceleration measurements. The lower horizontal line is the 
limit on .Mbh imposed by the absence of a rise of the stellar density 
expected inside the sphere of influence of the black hole (Ferraro et 
al. 2003a). 

As regard to the case of the single intermediate-mass 
black hole .Mbh in NGC 6752, we impose a mass of 
~ 500M Q . This value is derived according to the re- 
lation log(A4 B H/M Q ) = a + /31og(cr/200kms- 1 ) with 
a = 8.13 ±0.06 and (3 = 4.02 ±0.32 (Tremaine et al. 2002) 
extrapolated down to a ss lOkrns -1 . The central value of 
the line-of-sight dispersion velocity a in NGC 6752 has 
been estimated to be between 2.1 — 9.7 km s -1 (Dubath, 
Meylan, & Mayor 1997). Recent Fabry-Perot spectroscopy 
of single stars in NGC 6752 has shown a flat profile with 
typical dispersion of ~ 7kms~ 1 within the central 1' (Xie 
et al. 2002), while proper motion measurements of stars in 
the central part of the cluster suggest a much higher value 
of - 9 - lSkms" 1 (Drukier et al. 2003). In Figure 2 we 
have drawn the .Mbh versus a relation and indicated the 
upper limit on the central under-luminous mass implied 
by the pulsar acceleration measurements, together with 
the upper limit on the black hole mass .Mbh < 1, 000M© 
imposed by the lack of a rise in the stellar density inside 
the core resolved down to a scale of 0.08 pc (Ferraro et al. 
2003a). 

3. BINARY-BINARY ENCOUNTERS 

3.1. The binary pulsar 

Here, we study the 4-body dynamics of two binaries 
interacting under Newtonian gravity. Four cases are ex- 
plored resulting from the combination of "projectiles" oc- 
curring in two flavors, and "targets" of two types. The 



projectile is labelled as [PSR-A,CO] when it coincides with 

3Mq is the most general upper limit on the mass of a neutron star that we considered as lower limit for the mass of the lighter black hole. 
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the observed binary pulsar (in its post-recycling phase). It 
is labelled as [NS,CO] when the neutron star, not recycled 
yet, is orbiting around a more massive main sequence star 
(in its pre-recycling phase). 

The binary pulsar [PSR-A,CO] is described as in 
D'Amico et al. (2002). PSR-A has likely experienced a 
phase of recycling and of orbital circularization (Bhat- 
tacharya & van den Heuvel 1991) that has driven the 
neutron star to spin at the observed period of 3.27 ms 
(CPG). PSR-A orbits around a companion star of mass 
m co = 0.2 M Q (a value derived adopting a mass of 1.4M Q 
for the pulsar and an orbit inclination in the range of 
60-90 degrees). The binary separation and eccentric- 
ity are a PS R-A = 0.0223 AU (P or b,PSR-A = 0.86 days) 
epsR-A < 10~ 5 , respectively (D'Amico et al. 2002). 

In the pre-recycling phase, the neutron star (NS) has a 
mass of 1.4M© and the companion a mass m co of 0.8M Q 
consistent with evolutionary scenarios of low-mass binary 
pulsars (Tauris & Savonije 1999). The initial semi-major 
axis of the binary [NS,CO] is of 0.03 AU corresponding to 
an orbital period P rb,B = 1-3 days, below the bifurcation 
point (CPG). The initial eccentricity is 0.7. Here on, the 
total mass of the binary will be denoted as tub (1.6 M & 
for [PSR-A,CO] , 2.2 for [NS,CO] ). 

3.2. The [BE, BE] or [BE,*] binary as target 

The target is either a binary composed of two black 
holes [BH,BH] (in one case a black hole-neutron star bi- 
nary [BH,NS] ), or an intermediate-mass black hole with 
a star [BH,*] orbiting around. 

In the first hypothesis, the initial semi-major axis asH 
of all our black hole binaries [BH,BH] is of 1 AU, un- 
less specified otherwise, and the initial eccentricity is 
esH = 0.7. The black hole binary is always sufficiently 
hard to provide the right recoil speed V^psr-a to [PSR- 
A,CO] (or [NS,CO]). Since this speed is or the order of 
30 — 40kms _1 , the separation <jbh fulfills the inequal- 
ity: obh < 6^BH(MBH/lOM0)(4Okms _1 /Up : psR-A ) 2 AU 
(see eq. [3] and CPG for details; ^bh is the re- 
duced mass of the black hole binary). The bi- 
nary is also wide enough to have a coalescence time 
row longer than ~ 10 Gyrs which implies obh ^ 

0.4 [( MB H/lOM )(M T /lOOM ) 2 (T GW /lOGyr)] 1/4 AU, for 
an eccentricity of 0.7. This is a necessary condition since 
the characteristic time-scale for dynamical friction tdf to 
drive [PSR-A, CO] from the current position to the core 
has been estimated ~ 1 Gyr (see CPG; Sigurdsson 2003). 
This time-scale imposes a lifetime for the black hole bi- 
nary comparable to the age of the cluster itself, unless the 
binary formed recently. As an example, for a target binary 
with (50Mq, 1OM ), and eBH = 0.7, the suitable interval 
for a BH is (0.3,6) AU. 

In the second hypothesis, the target is the central 
intermediate-mass black hole of A^bh ~ 500 M sur- 
rounded by a swarm of bound stars belonging to the cusp, 
the region of influence of the black hole extending up to a 
distance tbh ~ GMbh/o- 2 ~ 0.02A / 1bh,50o/o' 2 pc, where 
cr io is the stellar central line-of-sight velocity dispersion 
in units of 10 km s _1 . (Note that even in the absence of 
a stellar cusp, a large central black hole A^bh can eas- 
ily capture a star or exchange with a binary star in the 
cluster core, likely a neutron star or a stellar-mass black 



hole.) We selected a cusp star (or companion star) of 1M Q 
tightly bound to the large black hole, moving on a Kep- 
lcrian orbit with semi-major axis obh,* of 1 AU, and ec- 
centricity 0.7 (-Porb.BH = 16 days). The cusp star (CS 
hereafter) is well inside the so called critical radius r cr i t 
where the flow of bound stars percolating across the loss 
cone peaks. At distances r < r cr i t relaxation leads to 
relatively small changes in the integrals of motion of the 
stars over a orbital period (Shapiro & Lightman 1977), and 
avoid prompt capture by the hole, the loss cone remaining 
empty. Thus, conditions q,bh,* < r cr jt and e<0.7 for a CS 
guarantees stability of the orbit during th e charac teristic 
time of the dynamical interaction T cnc ~ \J aBn,*b / a with 
[PSR-A, CO] (or [NS,CO] ), where b is the impact parame- 
ter of the encounter. In general, the following inequalities 
hold: P rb, psr-a < forb.BH ^ T eac < T rc i where r rol is the 
relaxation time inside the cusp. 



3.3. End-states 

Sampling of the initial conditions is carried on using the 
prescriptions outlined in Hut & Bahcall (1983), and Sig- 
urdsson & Phinney (1993). 

1 | 1 , , , 1 , , , 1 , n 



- [PSR-A, CO] out [PSR-A.CO] in - 




X (A.U.) 



Fig. 3. — [PSR-A, CO] is impinging onto the [BH,BH] binary of 
(50Mq,10Mq) and acquires the thrust to reach the halo of NGC 
6752. 

The integration method, a fourth order Runge Kutta 
scheme with adaptive step size and quality control, main- 
tains the accuracy on total energy and angular momen- 
tum conservation up to AE/Eq ~ 10~ 9 — 10~ 8 and 
AJ/J ~ 10~ n - 10~ 10 , respectively. After a close 
dynamical interaction, typical values of the fractional 
change of the binding energy of the black hole binary 
(APbh/Po.bh) are of the order of - 10~ 2 - 10~ 3 (here 
Eq,bh = G-Mbh to bh/2<2bh)- During integration, a 4-body 
to a 2-body switch is performed in the asymptotic regions, 
both for the incoming and outgoing states. The code has 
been tested against Sigurdsson & Phinney (1993) in the 
case of 3-body scattering. A trial 4-body experiment for 
a reference case (the 50M Q , 10M Q binary) has been com- 
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pared with the output of STARLAB 3 thanks to the coop- 
eration of Alessia Gualandris and Simon Portegies Zwart. 
We found agreement between the two outputs within the 
statistical uncertainties involved in the simulations 4 . 

Impact parameters and relative velocities of the two in- 
teracting binaries have benn selected to ensure an effec- 
tive energy exchange between the two binaries leading to 
post-encounter states that can be classified into five main 
groups: pure fliesby (FBs), ionizations (IONs), quartets 
(Qs), unresolved encounters (UEs) and exchanges (EXs) 
(the last two being very rare). Ionizations can further be 
sampled in (i) resonant ionizations in which [PSR-A, CO] 
is dissociated (their relative separation >20apsrt-A), but 
the two stars are bound to [BH,BH] (or [BH,*] ) (within a 
separation relative to the center of mass of black hole bi- 
nary smaller than 30cibh); (ii) ionizations with the escape 
of the two stars; (iii) ionizations with formation of a hier- 
archical triple in which PSR-A (or CO) escapes to infinity 
leaving CO (or PSR-A) bound to the [BH,BH] (or [BH,*] ) 
binary. These triplets may end either in stable or unstable 
states (Hut 1993; Mardling & Aarseth 2001; we return to 
these systems in §6). Quartets (Qs) correspond to cases 
where the total energy relative to the centers-of-mass of 
the binaries is negative 5 . UEs are very tight Qs occurring 
when the separation of [PSR-A, CO] relative to the black 
hole(s) never exceeds 30<2bh after 10 8 time steps. Integra- 
tion is aborted when Qs or UEs appear. On the contrary 
IONs and FBs are integrated until the outgoing star(s) 
reach the asymptotic state. As an illustration, Figure 3 
shows a FB for [PSR-A, CO] of interest to our studies. 

All end-states, i.e., FBs, IONs, UEs and Qs, are 
recorded to calculate relative probabilities, and for the 
case of FBs we produce post-encounter distributions for 
the recoil velocity, and eccentricity of [PSR-A, CO] (or 
[NS,CO] ). For [BH,BH] the distributions include IONs, 
UEs and Qs. All distributions are normalized to unity. 

3.4. Energy exchange and recoil velocities 

Similarly to the case of 3-body encounters (Hills 1983; 
Quinlan 1996), and given the large mismatch between the 
mass m B of the binary pulsar [PSR-A, CO] (or [NS,CO] ) 
and M T of [BH,BH] (or [BH,*]) we are led to quantify the 
mean energy exchange per scattering of the target binary 
([BH,BH] or [BH,*]) as 

A^bh r m B 

— sBH " 



focusing, is defined as 



E, 



0,BH 



(2) 



for the two 



where Air = A/bh + "Jbh or = A4 B h 
cases, and £bh is derived from our 4-body simulations. 
When the black hole binary hardens, most of its energy 
change will go into kinetic energy of the light projectile, 
giving a post-encounter recoil speed 6 

V 2 * G MBH ,„s 
^PSR-A ~ ?BH— (O) 
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the maximum impact parameter for 
energy exchange per scattering 
is estimated at start from 



-'oo.max 



the analytical expression of Sigurdsson & Phinney (1993), 
and later is determined more accurately from the numer- 
ical scattering experiments. Here, Voc denotes the pre- 
encounter relative velocity between the two centers of mass 
of the binaries; it is taken close to 7 — lOkms -1 , while 
V- 2 n = [2(mB+M T )(#o,BH + £o,B)/(m B M T )] is the veloc- 
ity necessary to dissolve the 4-body system (Bacon et al. 
1996). 
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4. RESULTS 
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Fig. 4. — Post-encounter recoil velocity distribution of [PSR- 
A,CO] including only FBs. [PSR-A,CO] scatters off the (50 M Q ,10 
Mo) black hole binary, and off [BH,*]. 

As an example, we first compare the end-states be- 
tween [PSR-A, CO] and the stellar-mass black hole binary 
[BH,BH] of (5OAf , 10Af s ) with those resulting from the 
interaction of [PSR-A, CO] off the black hole plus cusp star 
system [BH,*] of (5OOM Q ,1M ). We recall that the initial 
semi-major axis of the binaries hosting the black holes (or 
the hole and the cusp star) is of 1 AU corresponding to a 
binding energy of 4.5 x 10 48 erg (note that [BH,BH] and 
[BH,*] have equal energy in this case). [PSR-A, CO] at the 
observed separation of 0.0223 AU has a binding energy 



o B h 

where again /ljbh is the reduced mass of either [BH,BH] 
or [BH,*]. The cross section, enhanced by gravitational 

J http: //www. ids . ias . edu/ starlab/ 

4 Alessia Gualandris and Simon Portegies Zwart kinly accepted to run a comparison model, 
of the two codes will be presented elsewhere. 

5 Qs can be stable (Hut 1993) but are likely to be perturbed by flying by stars in the dense star cluster (Bacon et al. 1996). 
6 To avoid further labeling we will denote with Vpsr.— a the post-encounter velocity of the center of mass of the binary in the two cases 
[PSR-A,CO] , and [NS,CO] . 



A detailed comparison between the outcomes 
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of 10 47 erg. Fi gure 4 shows the post-encounter distribu- 
tion of VpsR-A obtained, for the two cases, collecting only 
FBs (the relevant end-states for the description of [PSR- 
A,CO] ). 

The distributions peak at V^psr-a ~ 30kms 1 (for 
[BH,BH] ) and 45kms~ 1 , (for [BH,*] ) just around the 
value of the recoil speed necessary to propel the binary pul- 
sar in the halo NGC 6752. The distributions are remark- 
ably asymmetric and we quantified their dispersion around 
the peak value in Table 1. Table 1 collects post-encounter 
physical quantities, while Table 2 collects data on rela- 
tive occurrence probabilities for the different end-states. 
In general, FBs are associated to less-vigorous/weak-recoil 
dynamical encounters compared to IONs. For [BH,BH] , 
FBs amount to 66% of all events, while only 11% of the 
outcomes end in IONs (see Table 2). For [BH,*] , IONs 
amount to 27% of the events against the 45% of FBs. The 
single massive black hole, having a stronger tidal field and 
a stronger gravitational focusing (the overall field being 
closer to a monopole) imposes a higher frequency of IONs 
ending in a narrower velocity distribution for the FBs. 

The post-encounter eccentricity of [PSR-A,CO] is re- 
markably different in the two cases. Figure 5 shows the 
two distributions. In the case of interaction of [PSR- 
A,CO] with a [BH,BH] binary, the end-state eccentric- 
ities peak around ~ 2 • 10~ 5 , and so [PSR- A, CO] pre- 
serves the very low initial eccentricity. The interactions 
off [BH,*] are instead quite damaging, yielding a signif- 
icant change in epsR-A- The difference can be ascribed, 
again, to gravitational focusing. In fact, in the interaction 
of [PSR-A,CO] off [BH,BH] gravitational stresses onto the 
incoming binary pulsar are reduced (due to the smaller 
total mass of the system) and softened (due to the dipolar 
nature of the interaction with the two BHs), thus causing 
less damage to the internal orbits of [PSR-A,CO] . 
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Fig. 5. — Post-encounter distribution of the eccentricity e of [PSR- 
A,CO] including FBs only. [PSR-A,CO] scatters off [BH,BH] of (50 
M Q ,10 Mq) and off [BH,*]. 

Table 1 surveys results from a series of run carried 
on considering stellar-mass and intermediate-mass binary 



black holes. The recoil speed of [PSR-A,CO] increases 
with increasing /xbh- But for the most massive binary, i.e., 
the (200M Q , IOMq) case, the rise of V^psr-a is faster 
than that implied by the linear scaling with /ibh (eq. [2]) 
suggesting a dependence of £bh on the total mass of the 
binary. For the (2OOM , 1OM ) binary, we find that the 
peak value of the recoil speed of [PSR-A,CO] largely ex- 
ceeds the escape speed. It varies with (Zbh as 

/i AU\ - 44 

V p ,psr-a = ^ p ,psr-a(1 AU) , (5) 

V a B H / 

so we can reduce the magnitude of the recoil considering 
FBs off a much wider [BH,BH] ; at a separation clbh of 7 
AU, [PSR-A,CO] would receive the right pull. Note that 
this intermediate-mass black hole binary does not damage 
significantly the eccentricity of [PSR-A,CO] despite the 
fact that the mass of the heaviest hole (200M©) is not so 
far from the mass A4bh of the single hypothetical black 
hole. 

The close encounters with the binary pulsar exert ran- 
dom impulses on the [BH,BH] binary. Linear momentum 
conservation then imposes a recoil velocity to the center- 
of-mass of the [BH,BH] binary. Figure 6 shows the skewed 
distribution of the recoil velocity for the (50M Q , 10M Q ) 
black hole binary scattering off the recycled [PSR-A,CO] . 
Table 1 collects characteristic values of the recoil speed 
computed averaging over all FBs and IONs. The rms re- 
coil velocity of the binary black holes is clearly in excess to 
its equipartition value (as expected for a point mass in a 
star background). It will be moderated by frictional drag 
and we defer to §7 for a discussion (see Merritt 2002 for 
a study of the Brownian motion of a massive binary near 
equipartition). 

As described by equation (2), the target black hole bi- 
nary transfers its gravitational binding energy to the bi- 
nary pulsar (mainly in the form of kinetic energy of the 
center-of-mass of the projectile binary). In the last row 
of Table 1 we collect the values of £bh obtained averag- 
ing over all encounters, and in bracket the values obtained 
selecting only FBs. The significantly lower value of £bh 
when inclusive of all events is related to the formation 
of triple systems: The formation of a hierarchical triple 
requires loss in binding energy of the target binary that 
widens to incorporate the third star. 

The black hole binary transfers its internal or- 
bital angular momentum Jo,bh to the orbit of the 
outgoing binary pulsar (that in this context can be 
treated as a point-mass, having negligible internal 
angular momentum). The initial Keplerian angu- 
lar momentum of the binary black holes, Jo,bh ~ 
10 55 (AfBH,2oo)(mBH,io)(M T:2 io)" 1/2 (aBH,iAu) 1/2 gcm 2 s- 1 
(subscripts refer to units consistent with a [200M Q , 10M Q ] 
binary) exceeds by a factor ~ 10 the or- 
bital angular momentum of the pulsar Jo b ~ 
10 54 TO B ,i.6(&oo/30AU)(K o /10kms- 1 )gcm 2 s~ 1 , 'where 
the impact parameter ~ (O.SobhGMt/V^) 1 / 2 and 
Vqo define the initial unperturbed hyperbolic orbit of the 
projectile. We find that the fractional angular momentum 
change of the black hole binary AJbh/Jo,bh amounts to 
~ 0.01 — 0.1 (see Table 1). This is comparable or less than 
Jo,b/ Jo,bh so the binary pulsar can only marginally mod- 
ify its post-encounter orbital angular momentum vector. 
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A more massive binary black hole (with mass >1000Af Q ) 
would be necessary to produce a sizeable change in the 
angular momentum vector of the binary pulsar. Stars 
of mass near turn-off ((m) ~ 0.5-M©) could instead be 
affected (see Miller 2003 for a discussion on the possible 
origin of rotation observed in the core of globular clusters 
and attributed to a binary black hole). 



r] — i — i — i — | — i — i — i — | — i — i — i — | — i — i — i — \ — : 

_ y (50, 10) M e _ 

; PRE-RECYCLING 



(50, 10) M s 
" POST-RECYCLING 
p 




V BH (km s-i) 

Fig. 6. — The post-encounter recoil velocity of the [BH,BH] binary 
(50 Mq, 10 Mq) in the pre- and post-recycling scenario. 

4.2. Pre-recy cling binary-binary encounters 

In this section we explore binary-binary encounters be- 
fore the neutron star is recycled, studying two cases: a 
(5OM , 1OM ) binary, and the single black hole [BH,*] . 

As shown in Table 3 the thrust received by [NS,CO] is 
within the correct range for expulsion of the binary (pro- 
genitor of PSR-A) in the halo of NGC 6752, in either cases 
([BH,BH] or [BH,*] ). Despite the large difference in mass 
between A^bh and the mass Mbh of the heavier black hole 
in the binary, the peak velocities Vp^sR-A are similar since 
recoil is sensitive to the reduced mass of the system. 

FBs off the [BH,*] binary lead to a significant increase 
of the post-encounter eccentricity and to a change in the 
binding energy of [NS,CO] . We have noticed that 20% of 
these FBs produce a [NS,CO] binary with orbital separa- 
tion such that the companion star fills its Roche lobe. This 
hints in favor of a "collisionally induced" recycling process, 
as suggested by Sigurdsson (2003). (Note that recycling 
lasts for a time less or comparable to the lifetime of the 
binary pulsar in the cluster halo according to dynamical 
friction, making this picture consistent.) Orbital changes 
of [NS,CO] are much less significant for the [BH,BH] case, 
and "induced" recycling has a frequency of less than 10%. 

In the pre-recycling case, the higher frequency of IONs 
found causes a remarkable increase in the recoil velocity of 
the [BH,BH] binary that can reach peak values ~ 4kms~ 1 
(Table 3). As illustrated in Figure 6 for the (50M Q , 10M©) 
binary, the distribution of the recoil velocity displays two 
humps, the first associated to FBs, and the second to 
IONs. IONs can thus imprint significant non-equilibrium 
velocities on the black hole binary. 



5. THE STATISTICS OF THE ENCOUNTERS 

In the post-recycling scenario, FBs are statistically more 
frequent. In general, with increasing mass of the black 
hole binary and decreasing mass ratio tobh/A/bh, the fre- 
quency of IONs increases; from 8%-12% if the target is a 
[BH,BH] binary, up to ~ 30% in the [BH,*] case. Triple 
systems have a high probability of formation, among IONs, 
and they will be described in §6. Qs account for ~ 20—30% 
of the end-states; many of these systems are wide and 
loose so that the interaction with the gravitational poten- 
tial of the cluster causes their fission into two separate 
binaries. In the pre-recycling case, the statistics reverse, 
causing IONs to become more important than FBs with a 
frequency as high as 85%, for the [BH,*] case. FBs com- 
prise only 7% of the events in the black hole-cusp star 
hypothesis, and 23% in the [BH,BH] case. 

Having computed the frequency of FBs relative to all 
possible end-states, we can estimate the collision rate of 
the binary off the [BH,BH] or [BH,*] 



T^coii = ~ £ B H,FBs27raBHGM T -i7 — > ( 6 ) 

T coll Voo 

where tins is the number density of neutron stars imping- 
ing onto the target black hole(s) and / the fraction in bi- 
naries with a star. The density of neutron stars is largely 
unknown given the uncertainties in their retention fraction 
at birth, but its value is expected to be larger than the den- 
sity n of stars (Sigurdsson & Hernquist 1993). The value 
of / is even more uncertain. To circumvent the problem 
of estimating tins and / we list, in Table 5, the minimum 
density of recycled neutron stars / tins necessary to have 
an encounter probability comparable to that of detecting 
[PSR-A, CO] at the current position. Since the estimated 
time of dynamical friction tdf is ~ lGyr, this minimum 
density is computed imposing a rate TZ co n ^Gyr -1 and 
a black hole binary separation <2bh equal to the maxi- 
mum necessary to acquire a recoil velocity of the order 
of ~ 35kms . These minimum densities cluster around 
50 — 10 3 pc -3 , implying a ratio fn^/n which is less than 
that expected from dynamical arguments (Sigurdsson & 
Hernquist 1993). Note that in the case of pre-recycling, 
/tins can be much larger, since there are no constraints 
on the spin history of the neutron star. 

6. EXOTIC END-STATES 

A black hole binary in a cluster is a catalyst for the 
formation of exotic triple systems composed of the two 
black holes and a star. This is a consequence of the high 
frequency of ionization events in binary-binary encounters 
ending with the capture of one of the two binary stars. 
Table 2 and Table 4 give in bracket the fraction of triple 
systems that form over the total. These comprise from 
<5% up to 30% of all events in [BH,BH] binaries. The 
question is now whether the triplet that forms is long-term 
stable, or short-term stable, i.e., destined to lose one of its 
components after a few outer-orbital times (Mardling & 
Aarseth 2001). 

In our experiments we have focused attention to binary- 
binary encounters leading to the formation of triplets 
with an active pulsar (in the post-recycling case) orbit- 



8 



ing around one of the two black holes 7 . We then extract 
from the sample those triplets that fulfill the long-term 
stability criterion following Mardling & Aarseth (2001; 
their eq. [90]). The stable triplets are more frequently 
composed by the millisecond pulsar orbiting around the 
heavier black hole (forming the inner binary) with the 
lighter black hole (the outer binary) orbiting around the 
inner pair. The typical parameters of the inner binaries 
found in our simulations have orbital periods ranging be- 
tween 10 to 100 days corresponding to orbital separations 
of 200 to 10,000 light seconds. The orbital acceleration im- 
parted to the millisecond pulsar by the black hole would 
not hamper its detection by current deep globular cluster 
pulsar surveys (Possenti et al. 2003; Ransom 2003; Ja- 
coby 2003). Figure 7 shows the formation of a triplet for 
the (5OM , 1OM ) binary. Triplets with the millisecond 
pulsar and the light black hole as inner binary are long- 
term stable only if their orbit is very tight. The stronger 
acceleration suffered by the millisecond pulsar due to its 
proximity to the hole would induce very rapid changes in 
the apparent spin period along the orbital motion impliy- 
ing a strong bias against the discovery of such pulsars 6 . 

MSP-black-hole binaries have long been searched in 
globular clusters as they can form via 3-body exchange 
interactions between single black holes and binary pul- 
sars (Sigurdsson 2003). We here for-see the possibility of 
detecting a hierarchical triple MSP-black-hole-black-holc 
system that would provide a direct measure of the mass of 
the heaviest black hole. 
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Fig. 7. — [PSR-A,CO] is scattering off the [BH,BH] binary 
(50Mq, IOMq) and the ordinary star escapes. A triplet forms with 
the millisecond pulsar orbiting around the heavier black hole. The 
inner binary has a semi-major axis of 0.28 AU and eccentricity 0.97. 
In the outer binary, the light black hole orbits around the center-of- 
mass of the inner binary with a semi-major axis of 5.1 AU, a perias- 
tron separation of 1.53 AU and eccentricity 0.7. The triplet fulfills 
the long-term stability condition of Mardling & Aarseth (with a ra- 
tio of the periastrion separation of the outer binary to the critical 
separation given by eq. [90] equal to 1.2, so greater than unity). 

Considering the case of a single intermediate-mass black 



hole in the cluster core, binary-binary encounters of the 
type studied here can lead to the "capture" of a millisec- 
ond pulsar thanks to the interaction of the binary pulsar 
with the start belonging to the cusp. This is a possibility 
that will be studied in more detail, introducing a model 
for the cusp. 

7. LINKING PULSAR'S ACCELERATIONS TO THE 
INTERMEDIATE-MASS BINARY BLACK HOLE 
HYPOTHESIS 

The single intermediate- mass black hole TWbh ^500 can 
account for a portion of the unseen matter required for ex- 
plaining the spin derivatives of PSR-B, PSR-E (and per- 
haps of PSR-D). 

Alternatively, one may wonder whether a suitably lo- 
cated perturber such as a black hole binary can accelerate 
the two pulsars. (This would reduce the demand of a sig- 
nificant amount of under-luminous matter in the core of 
NGC 6752.) The perturber should imprint a line-of-sight 
acceleration GMt/1 2 ~ c\P/P\. Considering the value of 
\P/P\ = 9.6 ± 0.1 x 10~ 17 for the two pulsars (PSR-B 
and PSR-E), and a separation I ~ 0.03 pc (D'Amico et 
al. 2002), comparable to the projected distance between 
the two pulsars, we can infer the minimum mass that the 
[BH,BH] binary should have to impart the observed accel- 
eration; M T > 180M Q . Thus, the (200M Q , 10M g ) binary 
has the right total mass. Assuming an harmonic potential 
for the central region of the cluster with uniform stellar 
mass density (m) n of ~ 10 5 M pc -3 (where (m) is the 
mean stellar mass), the minimum recoil velocity for mov- 
ing the black hole binary from the center of the potential 
well to the pulsar projected locations (at r ~ 0.08 pc) is 

<^>--( i ^)( I5F M^) 1/ \ m .- 1 . 

(7) 

Table 1 and Table 3 show that a [BH,BH] binary of the 
required mass has typical Vbh ~ lkms" 1 (2kms _1 ) in 
the post- (pre-) recycling scenario, lower than the value 
required. However the black hole binary is subject to 
repeated close encounters with cluster stars before dy- 
namical friction drives it toward the center of gravity of 
the cluster (Merritt 2002). This depends on the com- 
parison between the dynamical friction timescale tdf ~ 
a -3 / [1.89G 2 (m)nMx] and the collision time off stars r co u ~ 
a /[2itGM^na-Bn]. Their ratio 




is above unity but uncomfortably close. Because of this 
non-equilibrium dynamics, an intermediate-mass black 
hole binary of (2OOM , 1OM ) may random walk up to 
the current position of PSR-B and PSR-E. This scenario 
is possible but remains fine tuned. 

8. BLACK HOLE BINARY EVOLUTION 

In this Section we address the issue of the survival of our 
stellar and intermediate-mass black hole binaries against 
hardening by stellar collisions. 



^Hierarchical triples composed of two black holes and a main sequence star or a pulsar will be studied in detail in a separate paper (Mapelli, 
Colpi, & Possenti 2003). 
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8.1. Binary hardening 

The last black hole binary that remains in the cluster 
core (Portegies Zwart & McMillan 2000) has a large cross 
section and interacts with the cluster stars. Thus, we ex- 
pect it to harden progressively ejecting stars and emitting 
gravitational waves; at the same time stars may be tidally 
captured. 
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Fig. 8. — Black hole mass growth, Mbh versus time t (years) in a 
cluster with background stars with dispersion velocity a = 7kms" 1 
and density n = 10 4 and 10 5 pc — 3 . The dispersion velocity of NSs 
and WDs is obtained, in conditions of thermal equilibrium, re-scaling 

the stellar dispersion velocity as <Tj = yj (m) /rrn a, where (m) is the 
mean stellar mass and i=NS, WD. All lines terminate at the time of 
coalescence of the two holes by emission of GWs. 

It is the heaviest black hole in the binary that preferen- 
tially eats (unbound) stars because of its higher gravita- 
tional focusing. This leads to a secular change in the mass 
and orbital parameters. We mimicked evolution, solving 
the equations for the binary separation obh and mass A/bh 
as a function of time 



da 



BH 



dt 



= -27r£ B HGa BH X! 



mini 64G 3 Mbh^bh-Mt 

— C 5 a 3(l - e 2)V2 ' 



dMi 



BH 



dt 



= 2nGM*Q m 



(9) 



(10) 



2/3 njri 
a 



where the sum over index i is inclusive of all stars, i.e., 
main sequence stars, white dwarfs and neutron stars, dis- 
tributed according to the Salpeter IMF, and in thermal 
equilibrium. The capture rate on Mbh refers to a cross sec- 
tion equal to 27rGA/BH?~i,t/of where r iit = ^(ATbhA*'!*) 1 ' 3 
is the tidal radius of the star considered. We assume that 
the entire star is eaten. Figure 8 shows Mbh(*) for two 
values of the stellar density n (treated as a constant in 
time). We found that a secular increase in the mass of 
the binary black hole becomes appreciable only if the ini- 
tial mass Mbh(0) > 20M Q , under the background con- 
ditions imposed. A (50Af Q , 10M Q ) or (30M Q , 10M©) bi- 
nary avoid coalescence and its separation remains <1 AU 



for an appreciably long time. An "intermediate-mass" 
mass black hole binary (such as that considered in §7 with 
A/bh^200A/ q at the start) should have already terminated 
its life or be on the verge of coalescing, ending in a single 
black hole, and, curiously, should have grown to a mass 
close to .Mbh. Clearly, binary evolution is rather sensi- 
tive to the background density which is itself evolving with 
time (a fact that we did not include in this simplified treat- 
ment). Thus, the dynamical evolution of an intermediate- 
mass black hole binary should thus be addressed using a 
full N-body code. Lighter "stellar-mass" black hole bina- 
ries are still in the hardening phase and the mass increase 
is not significant. If our binaries will coalesce by GW emis- 
sion, they will remain always far from the critical regime 
(corresponding to mass ratios tobh/ATbh around 0.385) 
at which the gravitational rocket effect comes into play 
(Fitchett & Detweiler 1984), and will avoid ejection from 
the cluster due to release of non-zero net linear momen- 
tum. 

8.2. Star ejection 

Since black hole hardening (occurring at a rate 
dasu/dt = —2n£ i BHG{m)na 2 /a) is followed by stellar ejec- 
tion, a correlation exists between the mass lost in stars 



Afdcflc 



M T ln / obh(0) 
£bh \ obh (t) 



(11) 



and the total mass Mt of the binary black hole (treated as 
a constant). In equation (11), cibh(O) refers to the maxi- 
mum semimajor axis at which the interactions are strong 
enough to eject stars. When cibh(^) ~ acw, A/deficit is the 
mass that needs to be ejected in order to drive the binary 
toward coalescence by GWs. If NGC 6752 still hosts a 
binary black hole, one may ask whether the "deficit" in 
stars created by their ejection from the cluster would alter 
appreciably the stellar density profile or velocity field in 
the central regions. 

If, similarly to what claimed in bright "core-ellipticals" , 
a binary black hole turns an otherwise power-law density 
profile into a core (or shallower lower-low) due to the ejec- 
tion of stars (see Ferrarese et al. 1994; Milosavljevic & 
Merritt 2001; Milosavljevic, Merritt, Rest, van den Bosch 
2002 for details), then one can verify if the " mass deficit" 
(defined as the mass in stars that would need to be re- 
moved from an initially power-law profile in order to pro- 
duce the observed core) correlates with the mass of the bi- 
nary black hole according to equation (11). De-projecting 
the power-law+core profile of Ferraro et al. (2003a; Figure 
6) we infer for NGC 6752 a deficit of 300 "missing" stars. 
This is close to the mass that would be ejected by one of 
our binary black holes (eq. [11]). Since the central relax- 
ation time is of the order of the ejection time, a flow of 
stars can refill an underlying power-law, but gravitational 
heating (Spitzer 1987) induced by star loss can counteract 
this process. Only a detailed cluster model with a core 
binary black hole can address this question in detail, and 
should be designed specifically for NGC 6752. 

9. CONCLUSIONS 

What have we learned from this analysis ? 
(a) A stellar-mass binary black hole of (50Af Q , 10Af©) 
can imprint the right thrust to PSR-A in a gentle flyby 



10 



leaving almost unperturbed its eccentricity. Fliesby off this 
binary can occur at an acceptable rate without imposing 
unrealistic conditions on the neutron star density inside 
the central 0.1 pc. This binary black hole long lives in the 
cluster. Lighter black hole binaries or binaries with a black 
hole and a neutron star have interaction rates smaller by 
a factor ~ 10 and can be excluded (see Table 5). 

(b) An intermediate-mass binary black hole of 
(200M Q , 10M Q ) is more aggressive on the binary pulsar, as 
it tends to impart large recoil speeds to any incoming par- 
ticle. Thus, it has to be sufficiently wide (obh ~ 7—10 AU) 
to propel [PSR-A,CO] (or[NS,CO] ) in the halo of NGC 
6752 at the speed requested. Such binary demands the 
lowest neutron star density to allow one ejection every bil- 
lion years (Table 5). The survival of such binary black hole 
in the cluster however is uncertain: The binary may have 
merged into a single black hole in its interaction with the 
background stars: its presence in the cluster implies cither 
recent formation, or hardening on a time scale compara- 
ble to the age of the cluster. This binary may become the 
seed upon which our hypothetical "single" intermediate- 
mass black hole has grown. 

(c) A single, intermediate-mass black hole of .Mbh 
500M Q is a possibility. The recoil speed of PSR-A falls in 
the correct interval and considering its weak dependence 
on Mbh (recoil is more sensitive to the reduced mass of 
the system) even a 1000M Q would fit in this picture. The 
gravitational encounter with PSR-A must occur before re- 
cycling as the binary hosting the neutron star suffers a con- 
siderable change in its eccentricity Circularization can be 
achieved later, during the phase of Roche lobe mass trans- 
fer: we find that the binary is left in an outgoing state 
such to favor "scattering induced recycling" . 

The considerations in (a,b,c) do not address the problem 
of the acceleration of PSR-B and PSR-E. Can we provide 
a self-consistent picture ? If there is a stellar-mass binary 
black hole in the cluster propelling PSR-A, this binary can 
not act as a local perturber on the two pulsars. It is too 
light. The "best" candidate for the local perturber hy- 
pothesis is our intermediate-mass binary black hole (the 
200Mq, 10M© binary). The binary is sufficiently light so 
it can still random walk across the core to produce the 
local acceleration of the two pulsars. If instead the accel- 
eration of PSR-E, PSR-E (and perhaps PSR-D) is caused 
by the overall effect of the cluster potential well, the sin- 
gle intermediate-mass black hole plays a role but some 
additional mass should be present in the form of collapsed 



remnants such as white dwarfs and neutron stars. 

Black hole binaries are special catalysts for the formation 
of "stable" triplets. The interaction of a binary pulsar with 
an intermediate-mass binary black hole can create an ex- 
traordinary system: a millisecond pulsar-black hole-black 
hole hierarchical triple (the companion star leaving the sys- 
tem) . The single intermediate-mass black hole may, on the 
other hand, host in its cusp a millisecond pulsar. In this 
perspective, timing measurements on this "planetary-like- 
pulsar" , either member of a triplet or of the cusp, would 
give the unique possibility of discriminating the nature, ei- 
ther single or binary, of the black hole(s), providing also a 
reliable estimate of their mass. This opportunity exists in 
globular clusters contrary to the case of the Galactic Cen- 
ter, where the signal from a millisecond pulsar belonging 
to the cusp of our central black hole would be completely 
smeared out by scattering in the interstellar medium. 

A binary black hole may steadily perturb its environ- 
ment (§ 8.2). The distinction between a stellar- mass and 
an intermediate-mass binary black hole is subtle and dif- 
ficult to separate out. Accurate evolutionary models of 
clusters hosting a binary black hole of various masses need 
to be developed as they can indicate the fingerprints left 
by such exotic binaries. They also may provide clues for 
discriminating between the single versus binary nature of 
the central black hole(s) that may inhabit NGC 6752. 

During the revision of this manuscript Ferraro et al. 
(2003b) reported the discovery of the optical identification 
of the companion star to PSR-A, a helium white dwarf. 
The derived cooling age (<~ 1.2 — 2.8 Gyr) suggests that 
the dynamical encounter responsible for the ejection of 
PSR-A into the cluster halo occurred after the neutron 
star was recycled or was triggered at the time of ejection. 
In the last case the binary should have suffered very rapid 
evolution. 
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Table 1 

End-states in the post- recycling scenario. 





(500 M , 


l.OMg) (50 M , 1.4M ) 


(30 M , 3M ) 


(1OM , 10 M s ) 


(5OM , 1OM ) 


(2OOM , 10 M ) 


"Vp.PSR-A ( kms ') 
b ( v PSR-A> (kms- 1 ) 


45+30 

55 -13 


27ttt 


35± 6 2 2 


18 + 25 
18 -10 
o, + 112 


30+ 10 

7 o + 7S 
'• 3 -42 


go +eo 

y -60 

101 + 116 


° V p,BH (kms" 1 ) 
d (V BH ) (kms" 1 ) 


2+°-3 
0.3t° 8 2 


5+ 12 


8+ 2 ' 5 

1 9+ 3 ' 5 
^•"-1.1 


4.9l 10 6 9 


O-S-0.5 
9 i+2.8 


7 +0 - 7 

u - ' -0.5 

0-8^0.3 


° V rms,BH (kms" 1 ) 


0.465 


1.51 


2.45 


6.83 


2.62 


1.01 



S.PSR-A (X10- 5 ) 5^- 10* 2tY llT ill 3+^ 5t 



8 <ff>PSR-A (X10- 


3) 40±« 8 


4 + 712 


„ + 175 
"-9 


. + 4082 


9 + 577 
Z -2 


-1+ 6 

i -716 


h (|f >BH (X10- 2 ) 


o-itS:S 




1 0+ 18 ' 4 




1.8±«;» 


7+12.5 
u ' ' -0.7 


'<$J>BH (X10- 2 ) 


-o.oili-.il 


g + 63 
u -27 


9 + 25 


9+16 
_Z -12 


9 + 9 
Z -8 


-0.3+. 2 ,; 2 


?BH (?BH, FBs) 


0.3 (2.4) 


0.1 (0.9) 


0.2 (0.8) 


0.4 (0.5) 


0.7 (1.1) 


0.9 (1.6) 



"Peak value of the end-state recoil velocity distribution of [PSR-A,CO] considering only FBs. b Average recoil velocity of [PSR- 
A,CO] considering only FBs. c Pcak value of the end-state recoil velocity distribution of [BH,BH] binary considering FBs and 
IONs. d Average recoil velocity of [BH,BH] binary considering FBs and IONs. °Root mean square recoil velocity of [BH,BH] binary 
considering FBs and IONs. f Peak value of the final eccentricity distribution of [PSR-A,CO], considering only FBs. s Relative 
average change of the [PSR-A,CO] binding energy normalized to its initial value (Eo), considering only the FBs. h Relative average 
change of the [BH,BH] binary binding energy normalized to its initial value (Eo) considering FBs, IONs and Qs. 'Relative average 
change of the modulus of the orbital angular momentum of [BH,BH] normalized to its initial value (Jo) considering FBs, IONs 
and Qs. '^gjj as denned in eq. [2] obtained averaging over FBs, IONs and Qs. £bH FBs obtained considering only the FBs. 
The end-state distributions of the physical quantities are highly asymmetric. To quantify the skewness we have introduced asymmetric 
values for their dispersions. Dispersions around the mean values are calculated considering those values which contain 34% of the 
total area in the left and right wings, respectively. Dispersion around the peak values are calculated considering those values which 
contain 50% of the total area descending from the peak. 

The 12,000 binary-binary encounters simulated have been subdivided as follows: 1000 (for the 5Q0Mq,1Mq binary); 1000 (50,1.4); 
1000 (30,3); 1000 (10,10); 1000 (50,10);7000 (200,10). 



Table 2 

Statistics of the outgoing states in the post-recycling scenario. 



(5OOM , l.OM ) (50M Q , 1.4M ) (30M Q ,3M Q ) (10 M Q , 10 M Q ) (5OM , 1OM ) (2OOM , 1OM ) 



FBs 0.452 0.623 0.645 0.867 0.657 0.594 

EXs 0.0 0.001 0.0 0.0 0.0 0.0 

IONs (Triplets) 0.273 (0.267) 0.052 (0.043) 0.083 (0.072) 0.059 (0.048) 0.107 (0.092) 0.121 (0.103) 

Qs s 0.272 0.313 0.247 0.072 0.233 0.277 

UEs 0.003 0.011 0.025 0.002 0.003 0.008 
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Table 3 

End-states in the pre-recycling scenario*. 





(500 Mq, 1.0 M Q ) 


(50 M©, 10 M ) 


v p,PSR-A ( kms *) 
\ v PSR-A/ V^ 111B i 


45+40 
62+ 43 

u -18 


95+90 
zo -20 

67+i° 4 

ul -34 


V p,BH (kms^ 1 ) 
(V BH ) (kms- 1 ) 
V rms,BH (kms" 1 ) 


4+ - 3 
— 4 

5+ 1 - 3 
U.O_ f) 2 

0.602 


3.811-0 

3.2±2;» 
3.51 


e p, PSR-A 


9+ ' 1 

u - y -0.6 


4+0-5 


<|f >psr-a (xio- 1 ) 


-1 6+ 3 ' 4 

1 -°-1.3 


d.o_ 4 7 


(|f >BH (xlO- 2 ) 


1+ 26 


-1+53 
1 -1 


(■j~)bh (xio- 2 ) 


— 01 + 2 - 37 


-1+13 
L -9 


?BH (?BH, FBs) 


0.2 (1.7) 


0.3 (1.0) 



* Tabulated values are as in Table 1. 
We have simulated on a number 4000 of binary-binary encounters 
for every type of binary black holes (Mbh,™bh m solar masses): 
3000 (for the 500,1 binary); 1000 (50,10). 



Table 4 

Statistics of the outgoing states in the pre-recycling scenario. 



(500M Q , l.OM ) (5OM , 10M q ) 



FBs 0.068 0.232 

EXs 0.0 0.0 

IONs (Triplets) 0.844 (0.669) 0.686 (0.341) 

Qs 0.087 0.079 

UEs 0.0003 0.003 



Table 5 

Minimum neutron star density* 



M, m (Mq) 


a max,BH ( A - U 


f n NS (#/ P c 3 ) 


500, 1.0 


1 


2.9-10 2 


50, 1.4 


1 


1.5-10 3 


30, 3 


2 


1.110 3 


10, 10 


3.6 


7.6-10 2 


50, 10 


6 


2.0-10 2 


200, 10 


7 


5.5-10 1 



* Minimum neutron star density /hns m the cluster 
core to ensure an interaction every billion years, i.e., 
over a time comparable to the lifetime of PSR-A in 
the halo. The values of /n^g refer to the pre-recycling 
scenario for a (5OOM0,1M0) binary and to the post- 
recycling scenario for all remaining binaries. 



